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Abstract: Three- loop vertex diagrams in HQET needed for sum rules for B -B mixing 
are considered. They depend on two residual energies. An algorithm of reduction of these 
diagrams to master integrals has been constructed. All master integrals are calculated 
exactly in d dimensions; their e expansions are also obtained. 
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1. Introduction 

The mass difference Am in B°-B° is determined in the Standard Model by the matrix 
element <B°\Q(fi)\B°> of the four-quark operator 

Q(f,) = j a r, r = i L1 a d L (i.i) 

(see, e.g., Q). This matrix element is traditionally written as 

<B°\Q(fi)\B°> = 2 (l + ^ <S°|J a |0><0|J Q | J B°> J B(^) , (1.2) 

where N c is the number of colours. Here the first part of the right-hand side is the value 
of the matrix element according to the naive factorization prescription (this part does not 
depend on //), and B(/i) describes violation of this prescription. The hadronic parame- 
ter B{n) can only be obtained by using some non-perturbative method, such as lattice 
simulations (see, e.g., [|J) or QCD sum rules ]|, |||. 

In the QCD sum rules approach, the correlator <jQj> is investigated, where j is 
a current with <£>°|j[0> ^ (axial or pseudoscalar). Contributions to the theoretical 
expression for this correlator can be subdivided into two groups: 

<jQj> =2(1 + ^) <jJ a ><J a j> + <jQj> ni . (1.3) 



The first term includes the leading perturbative contribution plus all corrections (pertur- 
bative, vacuum condensates) to the two two-point correlators <jJ a >, <J a j> separately. 
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It just gives the square of the sum rule for Only the second, non-factorizable part 
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contributes to the sum rule for B{n) — 1. Non-factorizable perturbative contributions first 
appear at three loops (one gluon is exchanged between the two two-point correlators). In 
general, their calculation is a very difficult three-loop problem with three energy scales 
(m?, p 2 , p\\ we suppose that q 2 = 0) which cannot be solved at present. Several terms of 
the expansion in p\ , p\ have been obtained [|j (this is a much easier single-scale problem) . 
There are also non-factorizable terms due to vacuum condensates. 

It is also possible to consider sum rules in the HQET framework (see, e.g., |], |7|). 
The QCD operators Q, j can be expressed via HQET operators; matching coefficients 
are calculable series in a s (mb). Correlators of HQET operators don't involve the scale 
mi,. Therefore, no large logarithms appear in perturbative corrections. On the other 
hand, derivation and analysis of HQET sum rules for 1/nib corrections is difficult (though 
not impossible). Calculations in HQET are technically easier. In particular, three-loop 
diagrams describing the leading perturbative contribution to the sum rules for B — 1 involve 
only two scales — two residual energies. Here we present the method for calculating such 
diagrams. Calculation of this perturbative contribution is very desirable, because it allows 
one to control the /U-dependence of B{u) — 1. 

2. Reduction 







Figure 2: Diagrams to which the topology c 
reduces 



Figure 1: Generic topologies 

Non-factorizable three-loop diagrams be- 
long to three topologies (Fig. Q). Four HQET /^\^y^~^\ 
denominators in Fig. ||c are linearly depen- 1 
dent; therefore, one heavy line can be killed, 
and this diagram reduces to those in Fig. g 
which are particular cases of Fig. [ijh 

Let the incoming and outgoing residual momenta be pi, 2 . The scalar integrals depend 
only on the residual energies uj\ 2 = Pi, 2 ■ v, where v is the heavy-quark velocity. In the 
case uj\ = u>2 they reduce to single-scale HQET integrals || (see also H |i~0|]). 

We need to consider two topologies. The first one is (Fig. |||) 

1 fUjNp d d hd d k 2 d d h 
I a (n i]mj]UJl ,u 2 ) = J , (2.1) 



Di = -2(k 1 -v + u; 1 ), D 2 = -2(k 2 ■ v + w 2 ) , D 3 = -kf , D A = -k 2 
D 5 = -kj, As = -(k 3 - h) 2 , D 7 = -(k 3 - k 2 ) 2 , 



2 ! 



N\ = —2k 3 ■ v , N 2 = -{h-k 2 ) 2 



where — iO is assumed in all denominators, rii and nij are integer, and rrij > 0. They can 
be reduced to master integrals using integration by parts [11]. A Mathematics program 
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Figure 3: Topology 1 



(R.N. Lee, unpublished, based on [pig] ) has succeeded in constructing an algorithm to reduce 
these scalar integrals to the following simple master integrals: 



Mi(wi,u; 2 ) 



M 2 



, M 1 (uj 2 ,lo 1 ) = 





and one difficult integral: 



M' 2 




M 3 





(2.2) 



(2.3) 
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Figure 4: Topology 2 
The second topology is (Fig. f|) 

i /• EL K ] d% d d k 2 d d k 3 



D x = -2(k! ■ v + o>i) , £> 2 = -2(fc 2 • v + oj 2 ) , D 3 = -2((ki + k 3 ) ■ v + ut) 
D 4 = -kf, D b = -k 2 2l D 6 = -k$, D 7 = -(fc 3 - k 2 f , 
iV 1 = -(fci-fe 3 ) 2 , iV 2 = -(A:i-fe 2 ) 2 . 



(2.4) 
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The same program has succeeded in constructing an algorithm to reduce these scalar inte- 
grals to the same simple master integrals fl2.2|) and one difficult integral 



Ma 




(2.5) 



3. Simple master integrals 

We consider the below-threshold region uj\ 2 < 0; expressions for other regions can be 
obtained by analytical continuation. The simplest master integrals are single-scale, or 
products of single-scale integrals: 



\3d-7 



= I 1 I 2 (-2u Jl ) 2d ~H-^ 2 ) d ~ 3 . 



where the n-loop HQET sunset is 

I n = T(2n + 1 - nd)T n (f - l) . 
Several master integrals reduce to the one-loop vertex with two residual energies 

U 3 

1 f d d k 



(3.1) 
(3.2) 

(3.3) 




Ui U\ j Tl 2 UJ 2 



I(ni,n 2 ,n 3 ; wi, w 2 ) 



Z7T 



*/2 J D^D^D^ ' 



(3.4) 



D 1 = -2(k-v + u 1 ), D 2 = -2{k ■ v + uj 2 ) , D 3 



It is p] 



n\, n\ + n 2 + 2n 3 - d 
ni + ™2 



I(n 1 ,n 2 ,n 3 ;u ll u 2 ) = I{ri\ + n 2 ,n 3 ) 2 F 1 i 

x (-2w 2 ) d " ni " n2_2n3 , 
where the HQET two-point integral is 

r(m + 2n 2 - d)Y (I - n 2 ) 



x 



I(m,n 2 ) 



and 



Naturally, 



r(m)r(n 2 ) 



UJ 2 



(3.5) 

(3.6) 
(3.7) 



I{n 1 ,n 2 ,n 3 ;uj 1 ,U2) = I(n 2 , n 1 , n 3 ; u 2 , ^i) , 
/(ni,n2,n 3 ;w,£«;) = I{n x + n 2 ,n 3 )(-2u;) d -^- n *- 2n * 
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Later we shall also need 



l( ni ,n 2 ,n 3 ;u,0) = I (n u n 2 ,n 3 )(-2uj) d ^-^- 2 ^ , 

r (| - n 3 ) r(d - n 2 - 2n 3 )r(ni + n 2 + 2n 3 - d) 



T( ni )T(n 3 )T(d-2n 3 ) 



(3.8) 



Several ways to derive ( |3.5| ) are discussed in |l 
Using this integral, we easily obtain 



Mi(ui,u} 2 ) = I 2 I(5 - 2d, 1, 1;cji,uj 2 ) = I 2 I(1, 5 - 2d, l;u 2 , Ui) , 
M 2 (u>i,lo 2 ) = if 1(3 - d, 3 - d, l;u)i, u 2 ) , 
M 2 (uj\,uj 2 ) = if 1(3 — d, 3 — d, 2; wi, CJ2) • 

4. Master integral M 4 

We were able to calculate a more general integral 



(3.9) 
(3.10) 
(3.11) 



J(ni,n 2 ,n3,n A ,n 5 ;uJi,uJ 2 ) 




(4.1) 



Substituting ( |3.5|) for the left one-loop vertex subdiagram, we have 



J(m + n 3 ,n 4 ) 



?7T 



d/2 



(-2wi) 



i— ni — ri3— 2ri4 



dfcn d^A? 



(_fe2)n 5 (_ 2 (fc 0+a ; 2 ))n2 



2-^1 



n\,n\ + ti-3 + 2n 4 — d 
ni + n 3 



Then we perform Wick rotation = ik,Eo and take the d d 1 k integral. The integrand has 
a cut from to +ioo; we deform the integration contour around this cut (&eo = i(~ Lj °2)z): 



I(ni + n 3 ,n 4 )r (77,5 



d-1 ' 

2 , 



cos 



7r l/2 2 d-2n 5 -i r ( n5 ) 

PCO dzz d - 2 " 54 „ /ni,ni+n 3 + 2n 4 -d 
x / — — — 2F1 



[n (f - n 5 )] (-2wi) d - ni -" 3 - 2n4 (-2cj2) d ~" 2 " 2 " 5 



xz 



/0 (z + l)™2 - - y nx+n 3 
This integral can be calculated in terms of two 3 F 2 functions, and we arrive at 

Tt s r(f-n 4 )r(f-n 5 ) 

J(ni, n 2 , n 3 , n 4 ,n 5 ;wi, u) 2 ) = r(n 4 )r(n 5 ) 

T(ni + n 3 + 2n 4 — d)r(n 2 + 2ns — d) 
r(n 2 )r(ni + n 3 ) 

n 3 , ni + n 3 + 2n 4 -d,d-2n 5 \ _j- n2 - 2ri5 



x 3 F 2 



n\ + n 3 , d - n 2 - 2n 5 + 1 



x x 
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r(d — n 2 — 2n§)T(n 2 + n 3 + 2ns — d)T(ni + ri2 + n 3 + 2n 4 + 2ns — 2d) 
T(n 3 )T(d - 2n 5 )r(n! + n 2 + n 3 + 2n 5 - d) 



x 3^2 



n 2 , n 2 + n 3 + 2n 5 - d, n\ + n 2 + n 3 + 2n 4 + 2n 5 - 2d 



n 2 + 2n 5 — d + 1, n\ + n 2 + n 3 + 2n 5 - d 

X (_2o; 1 ) 2c( ~ ni_ra2_ra3_2n4_2n5 



(4.2) 



Trivial cases are reproduced: 

J(ni,ra 2 ,0,n 4 ,n 5 ;u;i,u; 2 ) = ^( n l> n 4)I("2, «5)(-2wi) 



d—ni — 2n,4 



(-2wz) 



d— ri2— 2n,5 



J(ni,0,n 3 ,re 4 ,n 5 ;u;i,u; 2 ) = J(n 3 , n 5 )J(ni + n 3 + 2n 5 -d,n 4 )(-2u;i 



At = w 2 , the single-scale integral |8| 
in 0). 

Now it is easy to write down the master integral fl2.5|) 

M^{oji,uj 2 ) = hJ(l, 1,3 — d, 1, l;o>i,o;2) 



2<i— ni — n3 — 2ri4 — 2715 

is reproduced (its derivation is also discussed 



(4.3) 



Note that the first 3 .F 2 function in Eq. ( f4.2[) turns into 2 F\ when one substitutes n 2 = 1 in 
order to obtain M 4 . 



5. Master integral M 3 

This integral can be expressed as 

M 3 (^i, w 2 ) 

via the two-loop integral 



£1/(1,1,1,1,2- g;wi,W2) 



7(ni, n 2 , n 3 , n 4 , n 5 ; w 2 ) 




n 2 w 2 



(m'V 2 ) 2 



£3 



-2(fci • « + wi) , D 2 



-2(A; 2 • w + w 2 ) , 



-A; 2 , D A = -k 2 2 , D 5 = -(ki-k 2 f 
with non- integer 77,5, where 

r(n + 1 



n|) (f - 1) 



T((n + l)(f-l)) 



(5.1) 



I 1 ,(5-2) 



(5.3) 



is the n-loop massless sunset. 

In order to express M 3 in closed form, we can use the method of differential equa- 



tions [14, 15]. The differential equation for this master integral can be obtained by dif- 
ferentiating it with respect to iv\ and then applying the reduction rules obtained by the 
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Mathematica program. It reads 



dM 3 {uji,uj 2 ) 



3d - 10 



duj\ 

H(u)i,U) 2 ) = 



(5.4) 



2d-5 (3d -7) [(3d -8)^-^-14)^1 3d _ 7 

Ml(i02,Ul) —73 oX 2/ \9 J -3\.—*U2) 



2(oj 1 -u 2 ) 2 v ' 1 8(d-3)w|(wi -c; 2 ) 2 

2d -5 (3d-7)[(3d-8)u; 2 -(5d-14)a; 1 ] 3d _ 7 
— ^Mi wi.wj H — — — — -2 1 3 (-2uji) 



Using the explicit expressions for the simple master integrals, it is easy to check that 
singularities at uj\ = uj 2 cancel in H separately on the second and third lines in Eq. ( ]5.4| ). 
The general solution of this differential equation has the form 



M 3 (u x ,u 2 ) = M {ux) 



C+ f 1 dwM - 1 (w)fr(w,wj) 
J —00 



where 



Af (w) = {-2uf d ' 2 - 



term in 



is the solution of the homogeneous part of the equation (5.4). In order to fix the constant 
C, we consider the asymptotics of M 3 (uji,uj 2 ) when u\ — > —00 16]. Using the method of 
expansion by regions (see flTfl ), it is easy to determine that there is no 0{uJ^ d ^ 2 5 ) 
the asymptotics. Thus, C = 0, and we obtain 

M 3 (wi jW2 ) = 2{-2u 1 f d / 2 -\-2u 2 f d / 2 ~^ (f - 1) T(8 - 3d) f°° — 



x \ y 



4-3d/2 



2F1 



1,8 -3d 
6 -2d 



y 



3d/2-4 



" / 1,8 -3d 






2Fl { 6-2d 


1 - 


■s) 







1 8 ~ 3 ^ 1 N 



8 - 3d / 1 
~ 6 -2d V ~ y 



(5.5) 



Note that the rational terms in brackets are the two first terms of expansion of the corre- 
sponding 2 Fi with respect to its argument. Now, using the parametrization 



2F1 



1,8 - 3d 
6 - 2d 



l-t 



r(6 - 2d) 



3d)r(d 



6-2d v ' 



ds s r-3d (1 + s) 2d-5 



1 



1 



s(l-t) 



i + st i + s (i + sy 



(5.6) 



we can take the integrals first over y and then over s. Finally, we obtain 



M 3 (uji,uj 2 ) = 4(-2^) M - 10 r 3 (|-1) 



r(8-3d) ^ 3d _ 9 

2(d-3) 



3-^2 



l,d-2, |d-4 
|d-3,3d-8 



3r(9 - 3d) 
2(d-3)(3d- 10) 



3-^2 



1,10 - 3d, 5 - |d 
6 - |d,4 - d 



x 
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+ 



7iT(6 - 2d) 
(3d - 10)r(d - 2) sin(3vr(i) 

7iT(6 - 2d) 
(d-4)r(d-2)sin(vrd) : 



2^1 



d,7-2d 



Id 



d-3 771 

2-^1 



— 7 — 
" 2' ' 

6 2 



2d 



(5.7) 



It follows from the analyticity of Ms(ui,U2) in the region u\ 2 < that the above expression 
is analytical in the interval x € (0, +00). In particular, branching singularities at x = 1 
cancel. 

The integral M3 is a symmetric function of its arguments. This symmetry can be made 
explicit if we rewrite the integral over y of the terms in the last line of (|5.5|) as follows: 



dy 



l/.l: 



dy 



d(l/y)y 2 



and make the replacement y — > 1/y in the second integral. Then, using the same parametriza- 
tion ( |5.6[ ), we obtain 



M z {uj u u 2 ) = (-2u; 1 ) 3a! / 2 - 5 (-2w 2 ) M/2 - 5 r 3 (d/2-l) 

r(|d-4)r 2 ( 5 -|d)r(2-f) 



(d-3)r(d-2) 



+ 2- 



3d) 



; 4~3 rf /2 3i?2 



+ 



d-3 
47rr(6 - 2d)x 3d / 2 ~ 5 



l,d-2,§d- 
|d- 3,3d- 8 



(3d - 10)r(d - 2) sin(3vrd) ' M 6 

o r(8-3cQ 3d/2 _ 4 /l,d-2,|d-4 
d-3 32 ' 



|d- 3,3d- 8 



+ 



4vrr(6 - 2d)x 



5-3d/2 



(3d - 10)r(d - 2) sin(3vrd) 




■|d,7 
6 - 



(5.8) 



We have performed two crucial checks of the above expressions for M3. The first 
check is due to the fact that at uj\ = 002 the integral M3 reduces to the known single-scale 
integral [filf . Though our representations do not literally coincide with those in filf , we 
have been able to check the perfect numerical agreement. 






Figure 5: Regions: thick lines are hard (momenta ~ loi), thin lines are soft (momenta ~ LJ2) 



The asymptotics of M3 (|5.7| ) at x — > can be also obtained by using the method of 
regions |17] for 1(1, 1, 1, l,n;u;i,a;2) with n = 2 — d/2 ([O]). There are 3 regions shown in 
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Fig. H The region a gives the first term in ( |5.7j ) ; b — the fourth term; and c — the second 
and the third ones (this is clear from the powers of oji, 2 ) • 

In the region a, we expand 1/D\ ( |5.2| ) in k\ ■ v. Then we calculate the left massless 
loop (lines 3 and 5) with the numerator {k\ -v) 1 (see, e.g., [17], eqs. (A. 11), (A. 12)). In the 
numerator of the remaining HQET integral, powers of 2A;2 ■ v may be replaced by powers of 
— 2^2, because integrals in which the denominator D2 cancels are zero. We obtain a series 
in x whose coefficients are finite sums. We have checked that a few terms in this series 
agree with the expansion of the first term in (|5.7|) . 

In the region b, we expand I/-D5 fl5.2j) in k 2 . Then we calculate the left (hard) HQET 
loop with a numerator (see |§, eq. (2.13)), and finally the right (soft) HQET loop (it also 
has numerators). The coefficients of the resulting series are finite sums. We have checked 
that a few terms in this series agree with the expansion of the fourth term in J5.7|) . 

In the region c, we expand I/D2 (j5l^) in ll>2- 



-2wi) 



2d-2n-6 



1=0 



where 

I(n h n 2 ,n 3 , n 4 ,n 5 ;uj, 0) = I ( ni , n 2 , n 3) n 4 , n ^- 2 u;) 2d - n ^- 2 ^- 2n ^ 

Using integration by parts, we obtain 

/ (1, 1 + 1, 1, 0, n + 1) - /„(!, I + 1, 0, 1, n + 1) 



/ (l,i + l,l,l,n) 



n- 



n 



where 



Io(ni,n 2 , 0, n 4 , n 5 ) = I(m, n 5 )I (nx + 2n 5 - d, n 2 , n 4 ) 
J (ni,n2,n3,0,n 5 ) = I{n 2 ,n b )I Q {n l ,n 2 + 2n 5 - d,n 3 ) 



(see (|3.q )). The contribution of the region c is thus 
r (| - 1) T (| - n - 1) T(2n + 6 - 2d) 



(d-n-3)r(n) 
r(2d - 2n - 5)r(2n + 3 - d) 



1 



d-3 




■ 2F1 



n + 3 — d, 2n + 3 — d 
n + 4- d 



x 



l,n + 3-d,2n + 6-2d 
4-d, n + 4 — d 



Substituting n = 2 — d/2 and multiplying by G% (see (5.1)), we reproduce the second and 
the third terms in (|5.7p. 



6. Conclusion 



We have considered scalar loop integrals needed for the perturbative part of HQET sum 
rules for B — 1. The sum rules will be considered in a future publication. The width 



-9- 



difference Ar involves matrix elements of four-quark operators similar to fll.lj ) but with 
different Dirac structures. In higher orders in l/m&, similar operators involving derivatives 
appear. Matrix elements of such operators can also be estimated using HQET sum rules 
(operators with derivatives are very difficult for lattice simulations). 

More general classes of three-loop HQET vertex diagrams can be analyzed using the 
same method. Master integrals calculated here will be useful for such an analysis. 

We are grateful to A. A. Pivovarov for discussions of HQET sum rules for B°-B° 
mixing. 

A. Expansions in e 



We use the Mathematica package HypExp |L9| to expand the master integrals in e (d 
4 - 2e): 



Afi 



r 3 (l-e)r(l + 6e) 



3x(l — x) 



72e 2 (l - 2e)(l - 3e)(2 - 3e)(3 - 4e)(l - 6e) 

1 \36x(l - xf logx - 6 + 71a; - 141x 2 + 105x 3 - 27x 4 



2 

i(l-x) 



18x(l - xf (8L(x) - 4 log 2 x - 9 log x) - 4 + 63x - 78x 2 + 21x 3 
+ (1 - x) 9x(l - x) 2 (48Li 3 (l - x) + 16Li 3 (l - x' 1 ) - 41og 3 x + 36L(x) 



M 2 



- 18 log 2 x + 7 log x) + 2(2 - 54x + 69x 2 - 18x 3 ) 

(1 - 4e)r 3 (l - e)r 2 (l + 2e)r(l + 6e) 

36e 2 (l - 2e) 2 (l - 3e)(2 - 3e)(l - 6e)r(l + 4e) 

3 r 
— -(1 — x) 2 (1 — x)(l + x) log X + X 



e 3 + 



(-2wi) 



4-6e 



(1 - x) 2 (l + x 2 



+ 



8(1 - x) 3 (l + x)L(x) + (1 - 2x - 2x 3 + x 4 ) log 2 x 
— 2x(l — x)(l + x) logx + 8x(l — x) 

96x 3 (2 - x) Li 3 (l - x) - 96(1 - 2x) Li 3 (l - x" 1 ) 

+ 24(1 - 2x - 2x 3 + x 4 )L(x) logx + (1 - x) 3 (l + x) log 3 x 
+ 24x(l - x)(l + x)(L(x) + logx) 



e 3 + 



ML 



-3x(l-x) 2 (3 log 2 x + 32 
(1 - 4e)r 3 (l - e)r 2 (l + 2e)r(l + 6e) 



4-6e^.-3e 



2 6e 3 (l-2e)(l-3e)(l-6e)r(l+4e) 
1 



(-2wi) 4 " b£ x 



x + -(1 - xfe 



+ 



+ 



3xlog 2 x + 3(1 — x)(l + x) logx + 4(1 — x) 2 
48x (2 Li 3 (l - x) + 2 Li 3 (l - x" 1 ) - L(x) log x) 
+ 24(1 - x)(l + x)(L(x) + logx) + (1 - x) 2 (3 log 2 x + 32) 



(A.l) 



(A.2) 



- 10 - 



+ 



(-2w 1 ) 2 - 6e x- 3e , 
T 3 (l -e)r(l + 6e) 



(A.3) 



36e 3 (l - 2e) 2 (l - 3e)(2 - 3e)(l - 6e) 
1 



6x-3(l + 13x + x 2 )e 



2x (91og 2 x + 16vr 2 ) + 18(1 - s)(l + x) log x - 9(1 + xf e 2 

48x (12 Li 3 (l - x) + 12 Li 3 (l - x" 1 ) - 6L(x) log x - 28C(3) - 5vr 2 ) 
+ 9(1 - x)(l + x) (16L(x) - 31og x) - 9(1 + 13x + x 2 ) log 2 x 



Ma 



- 12(1 - 15x + x 2 ) 
T 3 (l -e)r(l + 6e) 



e 3 + 



(-2wi) 



2-6 £a ,-3e 



(A.4) 



— X 

2 



6x log x — 1 + 18x + x z 



24e 3 (l-2e) 3 (l-3e)(l-6e) 

f 6x 2 (3 log 2 x - 2vr 2 ) - 18x(l - 9x - x 2 ) logx + 2 - 45x + 96x 2 + 15x s 



+ 6 



2x 2 (24Li 3 (l - x) + 24Li 3 (l - x" 1 ) - 12L(x) logx - 4vr 2 logx - 60C(3) - 9vr 2 ) 
+ 12x(l - x)(l + x)L(x) - 3x(2 - 9x - 2x 2 ) log 2 x - 3x(5 - 18x - 5x 2 ) logx 



+ x(7 + 2x - x 2 



£ 3 + 



(-2wi) 



3-6e 



(A.5) 



where 



1 



L(x) = -L(x _1 ) = Li 2 (l - x) + - log 2 x . 



As it was mentioned above, all the master integrals are analytical in x G (0, +oo), and 
hence the coefficients in the expansions ( |A.1| )-( |A~5 ) are analytical, too. It is easy to see 
that M2, Mg, M 3 are symmetric with respect to u\ «-> L02- The series ( |A.4[ ), (A.5) at 
oj\ = u>2 coincide with the expansions of the single-scale integrals p8l |^, P . 
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